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By using a combination of heteroepitaxial growth, structure refinement based on synchrotron x-ray diffrac¬ 
tion and first-principles calculations, we show that the symmetry-protected Dirac line nodes in the topological 
semimetallic perovskite SrlrOs can be lifted simply by applying epitaxial constraints. In particular, the Dirac 
gap opens without breaking the Pbnm mirror symmetry. In virtue of a symmetry-breaking analysis, we demon¬ 
strate that the original symmetry protection is related to the n-glide operation, which can be selectively bro¬ 
ken by different heteroepitaxial structures. This symmetry protection renders the nodal line a nonsymmorphic 
Dirac semimetallic state. The results highlight the vital role of crystal symmetry in spin-orbit-coupled cor¬ 
related oxides and provide a foundation for experimental realization of topological insulators in iridate-based 
hetero structures. 


I. INTRODUCTION 

Since the discovery of topological states in semiconductors, 
research on materials with strong relativistic coupling of elec¬ 
tron spin and orbital momenta, the spin-orbit coupling (SOC), 
have been flourishing owing to the huge potential for such 
materials to support spin-based electronics m. One of the 
most important frontiers focuses on the integration of SOC 
into the rich many-body physics of correlated electron sys¬ 
tems, in the search for novel phases, for example, a topologi¬ 
cal Mott insulator, Weyl semimetal, quantum anomalous Hall 
effect, spin liquid, and unconventional superconductivity in 
transition metal oxides mi. 

As a model system, perovskite Srlr 03 plays a key role in 
both addressing the SOC-correlation interplay and searching 
for novel topological phases in correlated oxides. Specif¬ 
ically, while the SOC-enhanced Mott instability has been 
shown to induce an unusual insulating state in Sr 2 Ir 04 and 
Sr 3 lr 2 07 from the Ruddlesden-Popper series Sr n+ ilr n 03 n+ i 
dlU, the n = oo member SrIrC >3 remains metallic |9 - JJ |. 
This dimensionality-controlled insulator-to-metal transition 
was believed to be caused by increased bandwidth (8). But 
recent angle-resolved photoemission spectroscopy (ARPES) 
lfT2l indicated that the bandwidth of SrIrC >3 is even narrower 
than Sr 2 IrC >4 with a semimetallic Fermi surface. Band struc¬ 
ture calculations mm have showed that, although SOC 
does aid in reducing the density of states at the Fermi level, 
full gap opening is prevented by band crossings at a Dirac 
nodal ring in the U — R — X plane originating from the mir¬ 
ror reflection of the crystalline Pbnm symmetry, rendering 
a topological semimetallic state. Moreover, it was suggested 
that lifting the Dirac degeneracy by breaking the mirror sym¬ 


metry would lead to various topological surface states in dif¬ 
ferent heteroepitaxial superlattices Ennm Experimen¬ 
tally, an ARPES study on epitaxial thin film samples observed 
a gap between Dirac bands rather than a nodal ring rm . n is, 
however, unclear how the symmetry protection and electronic 
structure would respond to epitaxial strain, which is essential 
prior to the realization of hetero structures. In fact, because 
the perovskite polymorph is not the thermodynamically stable 
phase of SrIrC >3 and bulk synthesis yields only poly crystalline 
samples (MU, epitaxial thin films are featured in the investi¬ 
gations of SrIrC >3 f8l fl8l42ni . including the two ARPES stud¬ 
ies Emo mentioned above, highlighting the necessity for 
understanding the influence of epitaxy on the symmetry. 

Here, using heteroepitaxial growth, structure refinement 
based on synchrotron x-ray diffraction and first-principles cal¬ 
culations, we show that an epitaxial constraint imposed by 
the substrate breaks the symmetry that protects the Dirac line 
nodes, and opens a Dirac gap. In particular, the results re¬ 
veal that the Dirac degeneracy is lifted without breaking the 
mirror symmetry, which is contrary to the previous theoretical 
proposals and indicates alternative possibilities for stabilizing 
topological states. A symmetry comparison further shows that 
another symmetry operation, i.e. the n-glide plane, within the 
Pbnm space group plays a more fundamental role in protect¬ 
ing the Dirac line nodes in orthoperovskite iridates, placing 
the system’s topological band crossing in the category of a 
three-dimensional nonsymmorphic Dirac semimetal. 

II. EXPERIMENTAL DETAILS 

High-quality single-crystal thin films (^18 nm) of per¬ 
ovskite SrIrC >3 were grown by pulsed laser deposition [2211231 
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FIG. 1. X-ray reflectivity data (red) and simulation (black). Thick¬ 
ness fringes and the reflectivity signal extend to high-momentum 
transfers, indicating the low roughness of the film/substrate interface 
and surface. Simulations with the Parratt formalism [24] yield a film 
thickness of 18 nm and a rms roughness between 0.2 and 0.4 nm for 
the substrate/film interface as well as the surface. 

on (llO)-oriented GdScC >3 ( Pbnm ), which applies ^ 0.6% 
tensile strain. In-house and synchrotron-based x-ray diffrac¬ 
tion experiments at beamline 6 IDB of the Advanced Photon 
Source confirmed the highly epitaxial single-phase perovskite 
structure and the film thickness (Fig. |T]). X-ray diffraction ex¬ 
periments along the truncation rods under grazing-incidence 
geometry were performed at SpLine (BM25, European Syn¬ 
chrotron Radiation Facility) with a fixed incidence angle of 
1° degree and an x-ray photon energy of 15 keV. Using this 
geometry, we performed crystal truncation-rod-scans along 
the surface normal ([110] direction) for more than 70 Bragg 
peaks. 


III. RESULTS AND ANALYSIS 

Symmetry. It is helpful to first illustrate the lattice symme¬ 
try operations of perovskite SrIrC> 3 , which has the orthorhom¬ 
bic GdFe 03 -type structure owing to the a~a~c + octahedral 
rotation [25i the same as the GdScC >3 substrate. This creates 
a four-formula unit cell from a \[2 x y/2 x 2 cell multipli¬ 
cation of the cubic perovskite cell. The [110] growth direc¬ 
tion is thus along a pseudo-cubic [100] c direction. The Pbnm 
space group, shorthand for P 2i/6 2i/n 2i/m, has a 6 -glide 
plane perpendicular to a (a be plane reflection plus a transla¬ 
tion along 6 ), an n-glide plane perpendicular to b (an ac plane 
reflection plus a diagonal translation within the ac plane), and 
a mirror plane perpendicular to c, as shown in Fig. [2] The 
6 -glide and n-glide operations are nonsymmorphic symme¬ 
tries, rendering Pbnm a nonsymmorphic space group. Tight- 
binding calculations suggested a Dirac nodal ring around the 
U point in the X — U — R plane originating from the mirror 


symmetry ED. In the following, we will experimentally show 
that this symmetry is preserved but the other two are broken 
in the strained film. 

X-ray diffraction and monoclinicity. In order to determine 
the unit-cell parameters, symmetry, and atomic positions of 
the Srlr 03 thin film, >70 different Bragg peaks were mea¬ 
sured in surface x-ray diffraction geometry (fixed incidence 
angle) using truncation rod scans along the surface normal 
[110] direction. The film diffraction signal is always found to 
have the same in-plane momentum transfer as the substrate, 
demonstrating the fully strained epitaxial growth. We also 
found intensity from the film near the ( 221 ) position of the 
substrate (Fig[3j)), indicating the bulk-like a~a~c + octahedra 
rotation (in fact a~b~c + as shown below) and discarding the 
other two; the [ 001 ] and [ 110 ] axes of the film are locked to 
that of the substrate within the surface plane and not rotated 
90°. 

Although the condition of epitaxy itself does not fix any 
interaxial cell angles, we found that a and (3 are 90° within 
the error bars by analyzing the diffraction peak positions. On 
the other hand, due to the lattice mismatch, the tensile strain 
introduces cell elongation along the [ 110 ] direction and com¬ 
pression along the out-of-plane [110] direction (Fig.[3]c), lead¬ 
ing to a deviation of 7 away from 90° and lowering the film 
lattice from Pbnm presumably to a monoclinic subgroup. In¬ 
deed, this angular deviation is reflected in the relative peak 
positions within each L-plane of the reciprocal lattice. For 
instance, if one uses the (400) and (040) peak positions to 
calculate those of (410) and (140) by fixing 7 = 90°, the 
simulated positions would clearly deviate away from that in 
the experimental scans and vice versa (Fig. -g). To account 
for this difference, 7 has to be increased by ^ 0.4°, consis¬ 
tent with the cell elongation along the [110] direction. Since 
7 % 90°, both the 6 - and n-glide operations become incom¬ 
patible with the lattice and the most plausible space group is 
monoclinic P112i/m (no. 11), which in Glazer notation lf25l 
corresponds to the a~b~c + pattern. 

To confirm this symmetry breaking, we simulated the com¬ 
plete set of truncation rods simultaneously and determined 
the structure factors. The simulation is required since a pre¬ 
cise determination of the diffraction intensity and position 
of the SrIrC >3 layer is hindered by the interference from the 
diffracted wave of the substrate (Fig. [3]). This interference 
leads to a deformation of Laue fringes, especially when the 
film and substrate peaks are in close vicinity to each other 
(Fig#). Using the thickness determined from x-ray reflec¬ 
tivity, we performed dynamical diffraction theory calculations 
1244 in order to determine both the position and structure fac¬ 
tor of the layer structure. These simulations nicely repro¬ 
duce the measured diffractions curves, including the afore¬ 
mentioned distortions seen in Figs. [3] and [4] Using this proce¬ 
dure, we obtain accurate values of the structure factor for all 
measured curves which are used to determine the atomic po¬ 
sitions in the unit cell. Using the structure factors, we refined 
the atomic positions in the unit cell and listed their values in 
Table |l] The results indicate that P112i/ra is indeed the best 
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FIG. 2. Schematics for the three symmetry operations of the Pbnm perovskite AB O 3 structure. For simplicity, only the BOq octahedra are 
shown. The blue octahedra with the black box denote the original Pbnm unit cell. The gray octahedra denote the unit cell after the reflection 
operation where the reflection plane is highlighted in orange. For the 6 -glide and n-glide symmetries, the solid/transparent octahedra denote 
the equivalent sites under each symmetry, respectively. The glide directions are denoted by dashed purple arrows. Since the mirror operation 
generates the enantiomorphic pair, the octahedra adjacent to the reflection plane become blue-gray due to superposition of the original and 
mirrored unit cells. 
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FIG. 3. a) The scattered intensity along the truncation rods in over 70 reflections obtained under grazing incidence diffraction geometry and 
plotted against the reciprocal lattice units of the GdSc 03 substrate, b) Truncation rod scan at the (221) Bragg peak, (c) The strain-induced 
monoclinicity between two orthorhombic lattice from the side view. Black, red, and blue show the substrate, unstrained and strained films, 
respectively. The orthorhombicity is exaggerated to illustrate the effect. Green arrows denote tensile strain, (d)-(g) Truncation rod scans at 
the (400), (040), (410), and (140) Bragg peaks, respectively. The diffraction signal from the substrate and thin film including Laue thickness 
fringes is experimentally observed (red) and reproduced by simulations (black). The peak position of the film and substrate are indicated as 
red and blue dashed lines, respectively. The gray dashed lines indicate the peak positions of (410) and (140) calculated based on that of (400) 
and (040) by fixing 7 = 90°. 


fit space group with 7 = 90.367(7)°. 

To rule out fitting artifacts, we performed correlation anal¬ 
ysis among the fitting variables, which confirms that the mon¬ 
oclinicity is only associated with parameters within the ab 
plane. Specifically, the unit-cell shape parameters, i.e. the 
lattice parameters (a, b , c) and unit-cell angles (a, /?, 7 ), 
were varied in order to simultaneously fit all measured crys¬ 
tal truncation-rod-scans. The angles and were found to be 
90.01(1)° and 90.00(1)°, respectively. The deviation from 90° 
is on the order of the error bar and similar deviations were 
found when fitting the diffraction positions of the orthorhom¬ 
bic substrate indicating that our the deviations are limited by 
our experimental resolution. On the other hand Fig. [3] illus¬ 
trates clearly that, in order to describe the layer peak posi¬ 
tion of selected diffraction peaks, a lowering of the symmetry 


to a monoclinic space group is required. To exclude possi¬ 
ble correlations between fit parameters as the source for this 
monoclinicity, i.e. 7 A 90°, we show the covariance ma¬ 
trix of the fit in Fig. [5{a). We find only values smaller than 
0.25 in the covariance matrix, indicating only weak correla¬ 
tions between the parameters. However, the strongest corre¬ 
lation is just found for 7 and the unit-cell parameters, a and 
b. In Fig. 0b), we therefore show the relative change of the 
fit error as a contour plot for the parameters 7 and a. A clear 
minimum is found for the values of a = 5.6120(5) A, and = 
90.367(7)°. Under no circumstances can the structure of the 
film be orthorhombic, since this would cause an significant 
increase of the fit error by more than 1500%. 

Density Functional Calculations. Based on the resolved 
crystal structure and symmetry, we computed the band struc- 
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FIG. 4. Selected crystal truncation-rod-scans (red) compared with 
dynamical diffraction theory simulation (black), (a)-(d) The data and 
simulations of the Bragg peaks (10-1), (301), (222), and (120), re¬ 
spectively. 


TABLE I. Atomic coordinates of Srlr 03 films on (llO)-oriented 
GdScCL substrates at room temperature. Space group is P112i/ra 
(No. 11), a = 5.6120(5) A, b = 5.5865(4) A, c = 7.934(2) A, 
7 = 90.367(7)°. 


Atom 

Wyckoff 

X 

y 

z 

Sri 

2e 

0.0047(2) 

0.9746(11) 

0.25 

Sr2 

2e 

0.4956(2) 

0.4741(14) 

0.25 

Irl 

2c 

0 

0.5 

0 

Ir2 

2b 

0.5 

0 

0 

Ola 

2e 

0.9420(15) 

0.4999(23) 

0.25 

Olb 

2e 

0.5571(7) 

0.9999(23) 

0.25 

02a 

4 / 

0.7179(23) 

0.2841(7) 

0.9709(12) 

02b 

4 / 

0.7790(15) 

0.7854(10) 

0.9738(27) 


ture for our strained film by local-density approximation 
(LDA) calculations with SOC included. The core and va¬ 
lence electrons were treated using the electronic configura¬ 
tions 4s 2 4p 6 5s 2 (Sr), 5d 7 6s 2 (Ir), and 2s 2 2p 4 (O); a 3x3x3 
Monkhorst-Pack k -point mesh f26l ; and a 600 eV plane-wave 
cutoff. Before performing calculations for thin films under 
strain, we calculated the band structure for perovskite SrIrC >3 
in bulk (strain free). Figure [ 6 ] shows the results from the LDA 
calculations without and with SOC. It can be seen that, in 
the absence of SOC, the degenerate t 2g orbitals form multi¬ 
ple intercrossing bands, giving rise to a ferromagnetic metal¬ 
lic ground state with a large density of states (DOS) around 
the Fermi level. When including SOC splits these bands, the 
obtained band structure consists of the empty e^-like bands 
found between 1 and 5 eV, and the t 2g states split by SOC into 
the occupied so-called J eff =3/2 states from -0.5 to -2.5 eV and 
the half-filled J e ff=l/2 states (0.5 to -0.4 eV). The J e ff = 1/2 


a 
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a b c a P y 
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FIG. 5. (a) Covariance matrix for the unit cell parameter fit showing 
that all cross correlations are smaller than 0.25. (b) Relative change 
of the fit error as a function of the lattice parameter and monoclinic 
angle 7. 


states form four intersecting bands owing to the four-formula 
unit cell and band-folding. The overall band structure and dis¬ 
persion are consistent with the other reported calculations for 
the bulk mm . 

Strain effect. Calculation using the resolved crystal struc¬ 
ture and symmetry of our strained film showed a similar over¬ 
all band structure (Fig. [7^). However, as shown in Figs. [7^ 
and b, we also found an ^ 5 meV gap between the upper 
and lower Dirac cones around the V point (the D point un¬ 
der monoclinic convention), indicating lifted Dirac degener¬ 
acy. Moreover, the Dirac gap transforms the band dispersions 
near the original Dirac nodes from linear to quadratic, i.e. a 
larger particle mass. To ensure the entire nodal ring is gapped, 
we surveyed multiple low-symmetry k -lines near the U point 
and confirmed that no point node exists (Fig. |7]c). This gap 
was unexpected because the Dirac line nodes were suggested 
to be a robust feature originating from the mirror symmetry 
m , which we found to be the only symmetry operation left in 
the film P112i/ra structure from the strain-free Pbnm space 
group (in addition to inversion). Note that the resolved struc¬ 
ture of our film is the only experimental input to our calcula¬ 
tion. To further validate this result, we examined our calcula¬ 
tions on the bulk Pbnm structure (Fig. [6]) and confirmed that 
the linear band crossing on the Dirac nodal ring is indeed well 
preserved therein (Fig.JTJi). 

Symmetry analysis. To reconcile this contradiction and un¬ 
derstand the symmetry protection on the Dirac nodal ring, we 
compared the symmetry breaking in our film to that of the 
Sr 2 lrRh 06 superlattice structure previously proposed for lift¬ 
ing the Dirac degeneracy ca. Note that, in the proposed su¬ 
perlattice, the mirror symmetry is designed to be broken by 
alternating the perovskite P-sites along the Pbnm c axis with 
Ir and Rh. Since this P-site alternation is only along the c 
axis, the 5-glide operation within the ab plane is preserved. 
However, the n-glide which involves a c-axis translation is 
required be broken simultaneously with the mirror symme¬ 
try. Indeed, the mass term introduced to achieve the mirror 
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Experimental structure 



FIG. 6. Band structures from LDA calculations without (left) and 
with (right) SOC. FM and NM stand for ferromagnetic metal and 
nonmagnetic metal, respectively. 


symmetry breaking also anticommutes with the n-glide oper¬ 
ator m. The resulting space group of a superlattice of this 
kind would be P2i/bll. Compared to our strained P112i/m 
structure, it is clear that the common broken symmetry is the 
n-glide. Therefore, to account for the lifted Dirac degener¬ 
acy in both structures, one can infer that the Dirac line nodes 
around the U point are protected by the n-glide symmetry op¬ 
eration, in addition to inversion and time reversal. 

We further corroborated this conclusion by performing cal¬ 
culations on structures of the P12i/nl and P2i/611 space 
groups, where only the n-glide and the 6-glide are left, respec¬ 
tively. The results showed that the Dirac nodal ring is indeed 
fully preserved under P12i/nl (Fig. [7^), demonstrating that 
the mirror or the 6-glide is not essential in the protection. In¬ 
deed, the nodal ring is also gapped under P2i/611 (Fig. HI) 
except for a pair of Dirac point nodes along the U—R (D—Eo) 
line. This is consistent with the calculation in Ref. ED where 
the nodal ring becomes nodal points in the S^IrRhOe su¬ 
perlattice. The difference is that the symmetry-removal was 
achieved here by unit-cell distortion rather than superlattice 
construction. Nevertheless, this observation shows that, in the 
absence of the n-glide, the Dirac points on this high-symmetry 
k -line can be protected by the remaining 6-glide symmetry. 
This is also consistent with the case of P112i/m discussed 
above; when both n-glide and 6-glide are removed, no Dirac 
points remain, and the nodal ring is fully gapped. 


IV. DISCUSSION 


This result leads to several direct consequences and impli¬ 
cations for understanding orthoperovskite iridates. First of 
all, one must consider the epitaxy-symmetry interplay when 
studying thin film samples. Because strain is inevitable in epi¬ 
taxial films, the crystal symmetry might be lowered and the 
Dirac degeneracy might be lifted. It is conceivable that the 



FIG. 7. (a) Band structures calculated for the strained film struc¬ 

ture given in Table [T] Red and blue bands highlight the hole-like and 
electron-like J e ff = 1/2 bands, respectively. By convention, the P-, 
X -, and P-points of the orthorhombic Brillouin zone are equivalent 
to the D -, P-, and To-points of the monoclinic one. The electronic 
LDA+SOC band structures near the P-point (P-point) for (b) the 
P112\Jm structure (determined by the strained Srlr 03 structure in 
Table|Tj), (d) the Pbnm structure of the bulk, (e) the P12i/nl struc¬ 
ture obtained by removing the 6-glide and the mirror, and (f) the 
P2 \/611 structure obtained by removing the n-glide and the mirror. 
Red and blue highlight the hole-like and electron-like bands, respec¬ 
tively. (c) Band dispersions along three different &-lines (inset) in the 
X — U — R(B — D — E 0 ) plane. Only results on a quarter of the 
X — U — R plane are shown here because the other three quarters 
are symmetrically equivalent. 


Dirac gap size will depend on the degree of deviation from 
the n-glide symmetry and increase with strain. Examples of 
strained SrIrC >3 films of this kind include that used in recent 
ARPES studies reported by Y. F. Nie et al. \\2\ and Liu et al. 
E). The former revealed the semimetallic bands with Dirac- 
cone like electron pockets, and the later reported a gap of 30- 
50 meV between the electron-like and hole-like cones near the 
U point with suppressed linearity, consistent with our conclu¬ 
sion. However, the epitaxy-induced symmetry-breaking was 
not considered in either study. A description of the band struc¬ 
ture taking into account the structure and symmetry is clearly 
requisite. 

Due to the protection by the n-glide operation, the Dirac 
semimetallic band crossing falls into the nonsymmorphic cat¬ 
egory of three-dimensional Dirac semimetals, which is now 
known as a fundamental class of topological semimetals 
127, 2j0. While other known nonsymmorphic Dirac semimet¬ 
als are characterized by point nodes l27l [29lh SrIrC >3 is dis¬ 
tinctive with its Dirac line node. Additionally, as discussed 
above, partial symmetry lowering may tune it into point nodes 
protected by another nonsymmorphic symmetry, the 6-glide 
operation. In fact, recently, we became aware of another the¬ 
oretical work by Fang et al. l30l . which also indicates that a 
four-band crossing nodal line on a k = n plane, such as that 
in Srlr 03 , can be protected by the nonsymmorphic symmetry 
in the presence of SOC. 

We can also further deduce the situations of other epitax¬ 
ial orientations, which we categorized here based on pseudo- 
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cubic orientations. Due to the nature of a glide plane, remov¬ 
ing it without complex superlattice structures requires strain to 
distort the reflection axis, i.e. changing a and/or 7 for the n- 
glide operation. The (100) c -orientation has three Pbnm vari¬ 
ants, (110), (lTO), and (001). It is obvious that (110) will in¬ 
duce the same monoclinic distortion and symmetry-breaking 
as (110). But the n-glide is preserved in the case of (001) be¬ 
cause [001] is collinear with the primitive axis and all angles 
may remain 90° except oblique surface lattices, e.g. LaAlC> 3 . 
The same is true for the (100)- and (OlO)-orientations, two of 
the (110) c variants. The third variant, the (lll)-orientation, 
is, however, opposite since the symmetry would become tri¬ 
clinic as strain is applied along all primitive directions. Fi¬ 
nally, (101) and (Oil) are the two (11 l) c variants relevant to 
many proposals of topological insulators in Srlr0 3 -based het¬ 
erostructures 10. For (101), since the 6-axis is within the 
surface and remains normal to the ac plane, the n-glide is 
preserved. It is, however, broken for (Oil) as this orthogo¬ 
nality is lost, i.e. a ^ 90°. The symmetry will decrease to 
P2i/611, gapping the nodal ring except for a pair of point 
nodes (Fig.JTjf). Therefore, the Dirac degeneracy removal sen¬ 
sitively depends on the epitaxial orientation. Whether Dirac 
gaps of different sizes can be open under the same amount of 
strain but different orientations would be an interesting sub¬ 
ject for further study. 

Finally, when the Dirac gap opens, the SOC-correlation in¬ 
terplay may significantly change as well. Note that Zeb and 
Kee d discovered from their calculations of bulk SrIrC >3 
that the larger SOC is, the weaker the correlation strength is, 
in sharp contrast to the insulating iridates, e.g. Sr 2 Ir 04 , where 
SOC is believed to assist correlation m This unusual be¬ 
havior in Srlr 03 is caused by the existence of the Dirac line 
nodes, which prevents full charge gap opening despite mod¬ 
erately strong correlation. However, the situation may change 
if the Dirac degeneracy is lifted; when the Dirac cones are 
free to renormalize under correlation, the SOC-correlation in¬ 
terplay may become synergistic. Therefore, while strain and 
film thickness may be used to engineer the bandwidth, elec¬ 
tronic, and transport properties, whether the n-glide symme¬ 
try and Dirac degeneracy are broken could lead to different 
responses. Recent resistivity measurements by several groups 
indeed reported strain- and thickness-dependent behavior and 
metal-insulator transitions in Srlr 03 films fTSHZI)- However, 
consensus on the underlying origin has not been reached, and 
the epitaxy-induced symmetry breaking remains to be com¬ 
pletely taken into account. 

VI. CONCLUSIONS 

In conclusion, by investigating the structural and electronic 
change in epitaxially strained Srlr 03 films, we showed that 
the epitaxy-induced symmetry breaking alone is sufficient to 
remove the Dirac degeneracy of the band crossing around 
the U point. We identified the symmetry operation protect¬ 
ing the Dirac degeneracy as the n-glide operation, which was 


hidden in the previous theoretical proposal for Srlr 03 -based 
topological insulators na. This finding places the system in 
the nonsymmorphic class of Dirac semimetals and provides 
new avenues to open the Dirac gap in much simpler strained- 
film structures and potentially stabilizing topologically insu¬ 
lating phases. Relieving this symmetry-protected band cross¬ 
ing may also reverse the SOC-correlation interplay. Our study 
thus highlights the critical role of epitaxy-induced symmetry¬ 
breaking, which must be taken into account in understanding 
and designing topological materials. 
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